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Abstract 

The neutrinoless double beta (Ois/3/3) decay is analyzed within the 
Minimal Supersymmetric Standard Model with explicit R-parity viola¬ 
tion (^pMSSM). We have found new supersymmetric contributions to 
this process and give the complete set of relevant Feynman diagrams. 
Operators describing 0 + —> 0 + nuclear transitions induced by the su¬ 
persymmetric interactions of the ^fpMSSM are derived. These operators 
can be used for calculating the 0 is/3/3 decay rate applying any specific 
nuclear model wave functions. 

The observation of neutrinoless double beta (Ou/3/3) decay would be a clear 
signal for physics beyond the standard model, since it violates lepton number 
by two units (for reviews see i. i). No definite observation of (V/?/? decay 
has been reported to date, but recent experimental progress has pushed the 
existing half-life limits of Ou/3/3 decay beyond 10 24 years and further progress 
can be expected in the near future 

Half-life limits on 0 u/3/3 decay are usually interpreted as limits on the effec¬ 
tive Majorana neutrino mass (see fig. 1). However, it is known for some time 
that there exist also other mechanisms which might induce 0 v(3f3 decays, like 
for example in left-right symmetric extensions of the standard model gauge 
group p|. In this paper we study contributions to 0 u/3/3 decay within super- 
symmetric (SUSY) theories with explicit R-parity breaking. R-parity (R p ) is 
a discrete, multiplicative symmetry defined as R p = (—\^ 3B + L + 2S ^ where S, B 
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and L are the spin, the baryon and the lepton quantum number. Conservation 
of R-parity in turn implies baryon number (B) and lepton number (L) con¬ 
servation. This symmetry has been imposed on the minimal supersymmetric 
standard model (MSSM) (for a review see |4|]). However, neither gauge invari¬ 
ance nor supersymmetry require R p conservation. The question whether or not 
R p is a good symmetry of the supersymmetric theory is a dynamical problem 
which might be related to more fundamental physics at the Planck scale. In 
general, R p can be either broken explicitly |J, 0 or spontaneously 0, || 
by the vacuum expectation value of the scalar superpartner of the _R p -odd isos¬ 
inglet lepton held || . Supersymmetric models with R p non-conservation (ft p ) 
have been extensively discussed in the literature not only because of their great 
theoretical interest, but also because they have interesting phenomenological 
1L 0) 0> and cosmological 0-|5| implications. One can expect that 


.^pSUSY models give a very natural framework for lepton number violating 
processes and particularly C0/3/3 decay [!2(]]. 

Let us start with a brief review of the model we use in the present paper. 

We follow the MSSM extended by inclusion of explicit .R-parity violating 
{]Jt,p ) terms into the superpotential . This model (^ P MSSM) has the MSSM 
held content and is completely specihed by the standard SU{3 ) x 577(2) x U{ 1) 
gauge couplings as well as by the low-energy superpotential and ’’soft” SUSY 
breaking terms J3J. The most general gauge invariant form of the superpoten¬ 
tial is |n|, [nj 

w = w Rp + w^ p . (1) 

The R p conserving part has the standard MSSM form 


W Rp = h L H x LE + h D HiQD - h v H 2 QU - pH x H 2 . (2) 

We use notations L, Q for lepton and quark doublet superhelds and E, U, D 
for lepton and up , down quark singlet superhelds; Hi and H- 2 are the Higgs 
doublet superhelds with a weak hypercharge Y = —1, +1, respectively. Sum¬ 
mation over generations is implied. For simplicity generation indices of helds 
and Yukawa coupling constants hjj, h D are suppressed. The mass-mixing 
parameter /i is a free parameter describing mixing between the Higgs bosons 
II i - II 2 as well as between higgsinos II\-II 2 . 

The Rp violating part of the superpotential ([!]) can be written as i, 0, 


— XijkLiLjEk + \ [j k LiQjDk + X"j k UiDjDk 


( 3 ) 


where indices i,j,k denote generations, and the helds have been defined so 
that the bilinear lepton number violating operators LiH 2 [6] have been rotated 
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away. The coupling constants A (A") are antisymmetric in the first (last) two 
indices. The first two terms lead to lepton number violation, while the last one 
violates baryon number conservation. Proton stability forbids simultaneous 
presence of lepton and baryon number violating terms in the superpotential 
[|0|] (unless the couplings are very small). Therefore, only A, A' or X" type 
interactions can be present. There may exist an underlying discrete symmetry 
in the theory which allows either the first or the second set of couplings. This 
discrete symmetry, called matter parity jlT], |T6| can be imposed ’ad hoc’ 


on the MSSM. It can be justified on a more fundamental level of Planck scale 
physics. Such a discrete symmetry is shown to be compatible with the ordinary 
SU( 5) H and ’’flipped” SU( 5) x 17(1) [Tjj] grand unification (GUT) scenarios, 
as well as with phenomenologically viable superstring theories ||18|| . 

Neutrinoless double beta decay, which is the main subject of the present 
paper, requires lepton number violating interactions. Therefore we bind our¬ 
selves to the lji p MSSM with lepton number violation (A ^ 0, X' ^ 0) and 
baryon number conservation (A" = 0). Apparently, 0 v/3/3 can probe only the 
first generation lepton number violating Yukawa coupling A' m because only 
the first generation fermions u, d, e are involved in this process. 

In addition to proton decay constraints on ^couplings there are also con¬ 
straints which follow from cosmological arguments requiring that the baryon 
asymmetry generated at the GUT scale is not washed out by B — L violat¬ 
ing interactions present in eq. (|3|). These cosmological constraints have been 
thought to affect all ^couplings A, A', A" -C 10~ 7 , making these models phe¬ 
nomenologically not interesting. These arguments, however, were proved to be 
strongly model dependent |T4|, JT5[] . Moreover it was shown recently that the 


cosmological bounds can be evaded in perfectly reasonable scenarios of matter 

M 


genesis 

The effect of ’’soft” supersymmetry breaking can be parametrized at the 
Fermi scale as a part of the scalar potential: 


Vsoft= £ rf\4>i\ 2 + h L A L H l LE + h D A D H l QD -huAuH 2 QU- (4) 

i=scalars 

—jiBHiH 2 + h.c. 


and a ’’soft” gaugino mass term 


1 r. 


Com = -- \M]_BB + M 2 W k W k + M 3 g a g c 


— h.c. 


(5) 


3 


























As usual, M 3) 2,1 are the masses of the SU(3) x £77(2) x f/(l) gauginos g, W, B 
and are the masses of scalar fields. Al, Ad, Ajj and B are trilinear and 
bilinear ’’soft” supersymmetry breaking parameters. Here, nig = M 3 stands 
for the gluino g mass. 

Now the model is completely specified and we can deduce the interaction 
terms of the ^yvlSSM - Lagrangian relevant for neutrinoless double beta decay. 

Let us write down these interaction terms explicitly. In the following we 
use the 4-component Dirac bispinor notation for fermion fields. 

The lepton number violating part of the Lagrangian can be obtained di¬ 
rectly from the W^ p superpotential part (|3]). It has the form 


- 


A' 


in 


(u L d 


R) 


+ (ul di) d r ■ 


°R 


da + (ez, v L ) d R 



( 6 ) 


-Vt 


+ h.c. 


To construct diagrams contributing to 0 v(3/3 decay we will also need R p con¬ 
serving gaugino-fermion-sfermion vertices. 

The Lagrangian terms corresponding to gluino L~ g and neutralinos C x in¬ 
teractions with fermions = {u,d,e}, q = {u,d} and their superpartners 
'ijj = {u, d, e}, q = { u , d} are |Q 


_ yO) 

C~ g = ~V2g^ ~ QrMr) + h.c., (7) 


C x = V2g 2 (tLiW'^LXii’L + CRii^RXii’R) + h.c. (8) 

2=1 

Here A^ are 3x3 Gell-Mann matrices (a = 1, Neutralino coupling 

constants are defined as [Q] 

CLi{i>) = -T 3 (ip)J\f i2 + tan^iy (T 3 (/0) - Q(^))Mn, (9) 

Cfli(V’) = Q(il>)taiL6 w ftJn. 

Here Q(il?) and T 3 (i/r) are the electric charge and weak isospin of the held 

V’- 

The coefficients N,j are elements of the orthogonal mixing matrix which 
diagonalizes the 4x4 neutralino mass matrix ||. The four neutralino mass 
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eigenstates Xi with masses m Xi 1 * have the field content 

X* = Mi b + M 2 iu 3 + Ms M° + M 4 ^2 0 . 


( 10 ) 


Recall again that we use notations IF 3 , B for neutral SU(2)l x U(l) gauginos 
and i? 2 , M for higgsinos which are the superpartners of the two neutral Higgs 
boson fields H 3 and H 3 . 

Having specified the Lagrangian interaction terms (BJ) - (§) one can con¬ 
struct diagrams describing the ^? p MSSM contribution to the 0 u/3/3 decay. The 
complete set of these diagrams is presented in fig. 2 . 

The supersymmetric mechanism of 0 u/3/3 decay was first proposed by Mo- 
hapatra ]20[ and later studied in more details by Vergados [^1|] . In these papers 
m , HU only three diagrams similar to those in Fig.2(a) were considered. In¬ 
stead of neutralinos Xii which are actual mass eigenstates in the MSSM, the 
consideration of refs. | 2 (J, |H] used Z-ino (Z ) and photino ( 7 ) fields in in¬ 
termediate states. Z and 7 can be mass eigenstates only at special values 
of parameters of the neutralino mass matrix. In general, these fields are not 
mass eigenstates. Using such fields in the intermediate states leads to neglect¬ 
ing diagrams with mixed intermediate states when, for instance, Z turns to 
7 due to the mixing proportional to the relevant entry of the neutralino mass 
matrix. The effect of mixing is taken into account completely in the diagrams 
displayed in fig. 2(a,b) with all neutralino mass eigenstates Xi involved. Note 
that diagrams in fig. 2 (b) were previously not considered in the literature. 

In the case of Ou/3/3 decay when momenta of external particles are much 
smaller then intermediate particle masses one can treat the interactions in Figs 
1,2 as point-like. A suitable formalism in this case is the effective Lagrangian 
approach. 

It is now straightforward to find the operators in the effective Lagrangian 
which correspond to the diagrams in Fig. 2(a,b). The result is 


£AL e = 2 , 


+ 


-eff 

e 


x 


= 8TTa 2 Xi U J2 


1 m Y , 

Xt L €-L 


e Li( e ) 


ml 


uld Ra )(uld R/3 )(e L e c R )+ (11) 


| 'P-(ulu^)(e L d^)(e L di) + 

TTi~ Til ~ 

U L dn 


1 The neutralino masses m Xi can be either positive or negative if the orthogonal matrix 

N is used in the diagonalization procedure [jf^. The sign of the mass coincides with the 

CP-parity of the corresponding neutralinos \i- 
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u L d R ) ( U 1 e< R) ( e L d R 


, ,j/3 > 


2 2 
777/j 777/~ /r 
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Here «2 = and a s = gf/(47r). 

The Lagrangian (pd]) has terms in a form which do not allow a direct 
application of the non-relativistic impulse approximation which is necessary 
for the standard calculation of the Op/3/3 reaction matrix element JI], [i2|. One 
should rearrange the right hand side of eq. (|TT]) in the form of a product of two 
colour-singlet quark currents and the leptonic current. It can be accomplished 
by a Fierz rearrangement procedure and subsequent extraction of color-singlets 
from the product of two colour-triplet and colour-antitriplet quark fields. The 
final result is 


r AL e =2 


0 ) 


m 


-i 


ri 2 

Ct p 

IT 

(v X e + v'g- V x f)JpsJps + VnJvaJvah (e(l + 75)e c ). 


(Vg + Vx)(JpsJps ~-WJtv) + 


( 12 ) 


The last term, corresponding to the heavy Majorana neutrino exchange di¬ 
agram in fig. 1, is included for completeness. The lepton number violating 
parameters are defined as follows 


Vg 


Vx 


V X e 


2na s A-f n m P 
9 G 2 F mA m~ q 

p drt y 
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m- 
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U L , 


7r«2 A 
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in 


6 G\m A 
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mp 


r • i m v . 
cLr 1 —1 


-Ri 


(■ d ) + e 


Li 


,U) 


= 27T cu 


a' 2 


111 


rrij 

ur 
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a >i 


777/, 


e L . 
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i =1 


Li\ 
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Xi 


Vg = 


Vxf 


47A'] 2 , j mp (rrip 


9 Gl- A 


F m d R m ~g 


m, 


TTO.2 An 


111 


m- 


UL / 

2 4 
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d-R 


, m 


UL , 


3 Glm\ 


m. 


e L . 


i=i m xi 




(13) 

(14) 

(15) 

(16) 
(17) 
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+ e Li (u)e Ri (d) ( ) + e Li (u)e Li (e) ' —^ 


m, 


U L . 


rriuL , 


W = 


ntp 


(18) 

< 771 N > 

where < > is the effective heavy Majorana neutrino mass (for definition 

see 0). 

Colour-singlet hadronic currents have the form 


Jps — w“(1 + 75 )d a , Jf 1 ' — u a a 7LV {l + 7 s)d Q , 

Jav = “Vfl - 7s)<t, 


(19) 

( 20 ) 


where a 711 ' = f[7 M ,7"]. 

Let us write down the 0 u/3/3 decay matrix element Tlovpp corresponding to 
the effective Lagrangian eq. 


Kovw = ^=m P 1 C 0 J {e(l + 7 5 )e c } x (21) 

(v-g + Vx) + { 7 he + v'g- v x i) ( F \nf\ J ) + VN( F \n N \iy . 

The heavy Majorana neutrino exchange contribution has been included in the 
last term of this equation. It corresponds to the last term of eq. (p~2|) . The 
normalization factor Cq u is defined below. 

We have introduced transition operators (for definitions see 0) which 
are useful for separating the particle physics part of the calculation from the 
nuclear physics one. The transition operators contain information about the 
underlying interactions at the quark level ( jT2| ) and quark states inside the 
nucleon. They are independent of the initial |/ > and the final < F\ nuclear 
states. To calculate the nuclear matrix elements in eq. one may use any 
nuclear model wave functions. In our case Qq,Qf and Qn describe transitions 
induced by quark colour singlet currents (|T9|)- (|20 ) in the first, second and third 
terms of eq. (]T2|). Diagrams in fig.1,2 with the intermediate states {W — N — 
W}, {u(d ) — x, g — u(d)} and {u(d) — x, 9 — d(u);e — X ~ e, q} contribute to 
operators Qn, Qq and Qj respectively. Calculation of the transition operators 
in the non-relativistic impulse approximation (NRIA) 0 requires the nucleon 
matrix elements of these currents. 

We derive nucleon matrix elements of the colour singlet quark currents 

H 


using results of ref. 


The necessary matrix elements are, 


< P(p)\ud\N(p') > = 


7 S) (9 2 ) ' N(p)T + N(p’), 


( 22 ) 
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Table 1: Nucleon form factor normalizations at q 2 


0 as calculated in [22 


Set 

b S 

^(3) 

r P 

t} 3) 

1 2 

rp( 3) 

1 3 

A) Bag Model 

0.48 

4.41 

1.38 

-3.30 

-0.62 

B) Non-relativistic 
quark model 

0.62 

4.65 

1.45 

-1.48 

-0.66 


< P(p)\u l5 d\N(p') > = Fp\q 2 ) ■ N(p)'y 5 T + N(p'), (23) 

<P(p)\ua^(l+ l5 )d\N(p l )> = N(p)(j^+ l -e^J pa y + N(p'), (24) 

< P(p)\uY(l ~ l 5 )d\N(p’) > = N(p)r (■ F v (q 2 ) - F A (q 2 ) l5 ) r + N(p'^5) 
where q = p — p' and the tensor structure is defined as 

= T (3) {q 2 )a pu + _JL_ W _ 7 V) + v _ a u P ^ (26) 

mp mp 


For all form factors Fv,A(q 2 ), Fs 3 p(q 2 )i T^\q 2 ) we take, following ref. 
dipole form 

^\{q 2 ) TfV) 


m 


FvAq 2 ) 


_ s,p\ 

fv,A F ( q 3) 


= 1 - 


g 

m 2 A/ 


-2 


,,, - | * 2, (27) 

s ,p(0) 7 ) 3 ( 0 ) 

with = 0.85GeV and fy ~ 1, f A ~ 1.261. Form factor normalization 
values Tg 3 p(0), t/ 3 \o) were calculated in ref. |22| and are given in Table 1. 
Using formulas (E2l)-(pB) we may derive the non-relativistic limit mp \p\ 


for nucleon matrix elements of the three combinations of the quark currents 
in (|12|). Keeping all terms up to order q 2 in the non-relativistic expansion we 
find the relevant transition operators Qq, 12^,12at for the case of two outgoing 
electrons in S-wave states, 


A 


-{a(?flp7v + o: a Agt,n + otyApi + o^^Igt 1 + (28) 

m P L J 


Vf = ttq(Ti = 0 ), 


m p f(fA 


u f 

n N 


m P 


Ja, 


UpTv — 12 


GT,N |; 


(29) 

(30) 


where partial transition operators are 


12 


GT,N ~ l^ T + T + a i 
*7U 


i?0 


— F n (xa), 


(31) 






























(32) 
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(*' 

T + 

> r (D 

T + 

*41 



E 

(i 

T + 

UD 

T + 

*41 



E 

(*' 

n 

UD 

T + 

*41 



E 

(*■ 

T+ 

> r (D 

r + 

*41 




Rn 


*? . 


F n {x a ), 

f Rq^ 


*J. 


Fi(x A ), 
— ) F 4 (x a ), 


R 


*7 . 


*7* 


a :j} 


Rn 


(33) 

(34) 

F 5 (x a )- (35) 


Here, Rq is the nuclear radius, introduced to make the matrix elements dimen¬ 
sionless (compensating factors have been absorbed into the phase space inte¬ 
grals HU). The following notations are used r t] = (T** — T^-), r %1 = \rij\, = 
r ij/ r ij, xa = rriArij. We also define the normalization constant in eq. (|2T|) 
Cq v = iir^j^m'X 2 . The above formulae have been written in the closure ap¬ 
proximation, which is well satisfied for large masses of intermediate particles 

i- 

The nucleon structure coefficients in (EH) are given by 


(Ay — 


Oly = 
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43)' 


/A / 

_ fm A \ 2 
\m P ) 

_ fm A \ 2 
\m P J 


(o) 
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(3)' 


/A 
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«A 
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(3)' 


T, 


(3) 




(3) 


T\ 


(3) 
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(o) 


2 T 


(3) 
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4 T 


(3)' 


i(3) 
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6 3 Tl (3) 3 jd 3) I 12 [f A 


exp = 


_ f m A V 
\m P J 
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( 0 ) 


IT, 


(3) 


2 T 


(3)’ 


+ -- 


12 3 t} 3) 3 Tj (3) 


1 

12 


n(3) 


Sa 


(36) 

(37) 

(38) 

(39) 


Here Fg 3 p = Fi? p ( 0), t/ 3) = t} 3 \ 0). Three different structure functions F t 
appear in eq. (|28|) (we use notations of ref. j2l|). They are given by 

(40) 


F n (x) = ^(3 + 3x + x 2 )e x , F 4 (x) = JH3 + 3a: - x 2 )e x , 


Fb(x) = %e 
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These functions are the analytic solutions of the relevant integrals over the 
intermediate particle momentum. They are analogons to the ’’neutrino poten¬ 
tials” for the case of light Majorana neutrino exchange |2|]. 

At this stage we point out that our formulas for the coefficients ( |36|) -( |39| ) 
of the transition operators (^8|)-(^9|) disagree with the corresponding formulas 
derived in ref. [EH]]. Particularly, in our case < 0 while in ref. [ET1] this 


coefficient is positive. This sign difference has an important consequence if one 
considers simultaneously the supersymmetric and the heavy Majorana neutrino 
contributions. Our formulas correspond to a constructive interference between 
these two mechanisms while formulas from ref. f2]| correspond to a destructive 
one. In the latter case both contributions can cancel each other and by a proper 
choice of A' nl in eqs. (|l^)-([T^) the matrix element of Oiyd/3-decay can be set 
to zero at any values of particle masses involved in the formulas. As a result 
the Quf3f3 -decay half-life limit would neither constrain supersymmetric particle 
masses nor that for Majorana neutrinos. Our formulas (P5|)-(PP|) always lead 
to certain constraints on these masses. Detailed discussion of this point will 
be given elsewhere. 

In conclusion we note that the results obtained in the present paper can be 
directly used for the extraction of constraints on the supersymmetric parameter 
space from the non-observation of Ozyd/3 decay. To do this one should employ 
nuclear wave functions calculated in a specific nuclear structure model. A 
subsequent paper will be devoted to this subject. 
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Figure Captions 


Fig.l Feynman graphs for the conventional mechanism of 0 vf3f3 decay by mas¬ 
sive Majorana neutrino exchange. 

Fig.2a Feynman graphs for supersymmetric contributions to 0 u/3/3 decay. 

Fig.2b New ” non-diagonal” Feynman graphs for supersymmetric contributions 
to 0 uf3(3 decay. 
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